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+^ Abstract 

S A formal proof to relate the concept of electromagnetic local density of states (LDOS) to the 
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I. INTRODUCTION 



The density of states (DOS), group velocity, and the distribution function are necessary 
for calculating various macroscopic quantities like specific heat, thermal conductivity, energy 
density, and radiation intensity. The local density of states (LDOS) is a generalization of 
DOS and, unlike DOS, is a position dependent quantity. Like the DOS, the LDOS depends 
on the type of carrier - electron, phonon, or photon. It is generally related to the Green's 
function of the appropriate governing equation (Schodinger equation for electrons, wave 
equation for long wavelength phonons, and Maxwell's equation for photons) and bound- 
ary conditions. In this paper, we are concerned with the photonic or electromagnetic LDOS. 

The electromagnetic LDOS is known to be related to the dyadic Green's function 
(DGF) of the vector Helmholtz equation [Ij. Quite often, the LDOS that is used in 
literature is related to the electric DGF and hence the electric field contribution [H 121 [3]. 
In free space, or vacuum, the electric field contribution and magnetic field contribution 
are equal in the absence of scatterers but that is not so in other circumstances. It was 
correctly pointed out that in addition to the electric field energy there is a contribution 
to the LDOS from the magnetic field energy and is related to the magnetic DGF The 
electric and magnetic DGFs are related to each other and will be discussed later in this 
paper. The reason for the usage of the electric DGF could be partially explained by the 
importance of the electric DGF in predicting the lifetime or decay rate of molecules in the 
vicinity of surfaces O E]. It is well known since the pioneering work of Purcell that the 
spontaneous emission rate of molecules is strongly affected by their vicinity and boundary 
conditions. A larger electric LDOS at the position of the molecule results in a shorter 
lifetime. Thermal near-field radiative transfer between a nanoparticle (or a dipole) and a 
large body, in the weak-coupling or first-order perturbation theory limit, can be explained 
in terms of the electric LDOS and hence the electric DGF [TJ |S]. Another related topic 
where the LDOS (both electric and magnetic) plays an important role is that of Casimir 
force between objects [HI ^U\. The Maxwell stress tensor in vacuum at thermal equilibrium 
can be expressed compactly in terms of the electric and magnetic DGF. 

The relation between the electromagnetic DGF and LDOS can be traced to Agar- 
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wal's work in which he used hnear response theory to express electromagnetic field 
correlations in terms of suitably defined response functions [IT]. These response functions 
are related to the electric and magnetic DGF. A different proof for the relation between 
LDOS and DGF using eigenfunction expansion of the DGF is given in pLj. The DGF is 
widely used in solving electromagnetic boundary value problems and there is a rich collection 
of works on this topic, of which only a few important ones are cited [121 [131 [IH [El US] . The 
singularity of the DGF is a topic of practical importance, especially for numerical solutions 
to scattering problems. Like the Green's function for the Laplace equation and the scalar 
wave equation, the DGF G{r, r') too exhibits a singularity as |r — r'\ 0. This singularity 
is stronger than that of the Laplace or scalar wave equation and the DGF behaves as 
\r — r'\^'^ as \r — r'\ 0. This singularity is generally dealt with by introducing a 
principal volume and the depolarization dyad, L, that depends on the shape of the principal 
volume. The DGF as well as the depolarization dyad are required to determine the electric 
and magnetic fields within a region containing sources and is topic of discussion in these 
works [m [T7|. Since the proofs for the relation between LDOS and the DGF or response 
functions in [H |TT] do not take this singular nature of the DGF explicitly into account, 
it is not clear whether they can be extended to relating LDOS to the DGF in material 
media, including lossy materials. In this paper, we define the electromagnetic electric and 
magnetic field energy density at any point, whether in free space or any dielectric material 
in thermal equilibrium at a temperature T. Temperature-induced thermal fiuctuations and 
quantum or zero-point fiuctuations of charges act as sources of electromagnetic radiation. 
The energy density has to be determined in a region containing the source and hence the 
necessity to take into consideration the L dyadic, in addition to the DGF. The expression 
for spectral energy density, U{r;uj,T) at angular frequency u can be separated into a part 
that is related to the average energy of a harmonic oscillator at temperature T, Q{uj,T), 
and the rest of which gives the LDOS, p(r; u) to give 



where 27rh is Planck's constant and is Boltzmann's constant. Though the L dyad is 
necessary to correctly define the field in the source region, the LDOS should be defined 
such that it does not depend on the L dyad since it should not depend on the shape of an 



U{r;u,T) = p{r;u)Q{u,T) 




(1) 
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arbitrary principal volume. 

The problem of electromagnetic energy density in material medium is a fascinat- 
ing topic and has a long history, especially with the recent development of electromagnetic 
meta-materials [iHl [191 EDI EH |22l [23] . Generally material media are described by frequency 
dependent electrical permittivity, e{uj), and magnetic permeability, /i(co'). This frequency 
dependence is what is known as temporal dispersion. One of the shortcomings of taking 
into account only temporal dispersion is that it cannot explain phenomena like natural 
optical activity or gyrotropy [24J. This was later explained by taking into account the 
weak dependence of permittivity on wavevector. This corresponds to taking into account 
length scales of the order of the molecules in the medium as opposed to a continuum 
theory, where there are no such length scales. Spatial dispersion of electrical permittivity 
or wavevector dependence of permittivity (when spatial dispersion is taken into account, 
permittivity and permeability are not independent quantities [23]) is important in investi- 
gating electromagnetic properties of plasmas and metals at low temperatures [201 [2S|- It 
is well know from Rytov's seminal work on thermal fluctuations that the energy density, 
as well as thermal radiation intensity, in absorbing media is infinite if spatial dispersion is 
not taken into account [191 . More recently, Tai and Collin have investigated the radiation 
from a Hertzian dipole immersed in a dissipative medium (no spatial dispersion) [27] . They 
found the total radiated power from the dipole to be infinite, which is what Rytov realized 
for electromagnetic radiation due to thermally fluctuating sources. Embedding the dipole 
in a cavity fllled with a lossless dielectric material ensures that the power radiated by the 
dipole is flnite but dependent on the size of the cavity. Generally the shape of the cavity is 
a sphere and the material of the cavity is free space or vacuum. The power radiated by a 
dipole at the center of the sphere is proportional (asymptotically) to R~^, where R is the 
radius of the cavity [271 [25] I2U] . One way of trying to overcome this problem is to include 
spatial dispersion in every material. Determining the DGF is complicated enough as it is 
with termporally dispersive materials, let alone ones with spatial dispersion. We will show 
in this paper that the LDOS at any point p{ro) can be split into two parts - one that does 
not depend on the location of and the other that depends on the position Vo'. 




(2) 
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Though po does not depend exphcitly on ro, it is an imphcit function through the 
permittivity and permeabihty at r^. For transparent materials, po is well defined and can 
be determined from the DGF for infinite, homogeneous space of the same material. For 
absorbing materials, that is not so. To ensure the finiteness of po the spatially dispersive 
nature of the material should also be taken into account. The second part, Psci^o), depends 
on scattering from boundaries, and can be determined by calculating the scattered part 
of the DGF, which exhibits no singularity as — ro| —>■ 0. Determining po for spatially 
dispersive materials will not be dealt with in this paper. As mentioned in Joulain's work 
[1], a proper definition for LDOS in material media, with possible losses is not known to 
the best of our knowledge. Our aim is to provide that. 



There exist multiple definitions in literature for energy density in a linear, absorb- 
ing dielectric [2n [23| 130] . The common used ones for electrical energy density are: 

{eU^, ii e = constant, ^(e) = 0, 

where U!: (u) is the total energy density and W: = (l/2)eo\E(r,u!)\'^ is the energy density 

I ot J 

due of the electric field, Q{z) refers to the imaginary part of z. A more general result is 
given by Loudon [21] and, later, Ruppin [23] for a dielectric material composed of Lorentzian 
oscillators . For a medium with electrical permittivity given by 

^^°'-)-^ + ^-^-,r. « 

the total electrical energy density is given by 

/ 2uif" \ 

uLM = [^ + <o)u; (5) 



where = e'^^ + ie'^^ . The expression for energy density in Eq. [5] coincides with Eq. [3] in 
the limit e'^^ — * 0. What we see from Eq. (j3]) and Eq. (jsj) is that while the expression for total 
energy depends on the type of the material, the field energy can be defined unambiguously 
and the total energy is a multiple of this field energy. Similar expressions exist for energy 
density of magnetic fields, [23J. We will show that the electric field energy density, f/J, 
and magnetic field energy density, f/™, can be written as t/J = pe {r,u) 6 {uj,T) and [/™ = 
Pm{r,uj) Q {uj,T), where Pe{r,uj) and pm{r,uj) are the electric and magnetic LDOS. The 



expressions for pe (f,^) and pm (f,^) thus obtained agree with the definitions of electrical 
and magnetic LDOS in vacuum as suggested in |3]. 

The paper is arranged as follows: We first derive expressions for electric and magnetic 
field in terms of the respective DGF in Section |TTj The boundary conditions for the DGF 
are obtained along the way. The fiuctuation-dissipation theorem gives us cross-spectral 
density of the thermal current density vectors. Though a significant portion of this section 
can be found in various sources, substantial detail is provided because of its importance 



in obtaining the eventual expressions for LDOS. In Section III we express the electric and 
magnetic field correlations in terms of the products of DGF. Finally, using Green's theorems 
for dyads, expressions for the cross-spectral densities of the electric and magnetic field 
vectors and expressions for the electric and magnetic LDOS are obtained. The vector- 
dyadic and dyadic-dyadic formulae used in this paper are given in Sec. |X] and Sec. |B} 



II. ELECTROMAGNETIC FIELDS AND DGF 

The configuration of the sources is shown in Fig. ([T|. The finite objects are all embed- 
ded in a host medium. Each of these objects, including the host are defined by relative 
permittivity, e {u; r) and relative permability, /i {u; r). Each object in itself is homogeneous 
and hence the permittivty and permeability only vary across a boundary between two ob- 
jects. Even though the host medium in general is vacuum, for the purposes of defining a 
temperature, and hence being able to use the fiuctuation-dissipation theorem, we take the 
material to have e {uj; r) and fi {uj; r) given hj 1 + i6 with 6 —>■ (the imaginary part could 
be different for e and /i but it does not matter finally). This idea of a tiny absorptive part 
that tends to is the mathematical version of the concept of carbon particle used by Planck 
in his treatise on heat radiation [311 page 44]. The spectral electric field and magnetic field 
are governed by the macroscopic Maxwell's equations: 

V X - iujfXof^H = - (6a) 

V xH + iuSosE = J"^ (6b) 

Eo and /io are the electrical and magnetic permeability of free space. and J"^ are the 
electic and magnetic current densities. The explicit r dependence of e and /i is suppressed 
henceforth. These equations can be converted to the vector Helmholtz equation for electric 
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and magnetic fields as: 

V xV X E -k'^E = iujjJof^J^ - V X J"^ (7a) 

V X V X - k'^H = V X J*^ + lUEoeJ'^ (7b) 

where k'^ = To invert these equations and express the electric and magnetic fields 

as integrals over the source regions, we make use of the DGF which also obeys the dyadic 
version of the vector Helholtz equation: 

V X V X ^,{r, r') - k^^e{r, r') = l5{r - r') (8) 

Here the V operates on r and k'^ is also defined at r. To obtain an integral expression for 



the electric field, we use Eq. (A4) with E{r) for F and Ge{r,r') for G. The volume of 
integration is V — Vs, where Vs is a small volume surrounding the point r' in order to avoid 
the singularity of the DGF, and V is the whole of space. To be able to express E as an 
integral over the sources, we need equality of the surface integrals of field quantities defined 
on either side of the boundary between two materials, which in this case is 



- n . {E{r) X V X G^{r, r')) + V x E{r) x Ge{r, r')) 



dS 



where h is normal to the boundary surface S. Using Eq. (Al) and Eq. (A2) the above 
expression is transformed to 



y (n X E{r)) . (^n X x V x Ge{r, '^'))) + (n x V x E{r)) . (^n x (^n x Ge{r, r'] 



dS 



Since the tangential electric {n x E{r)) and magnetic fields (n x H(r)) are continuous 
across the boundaries, the continuity of the surface integral demands the following boundary 
conditions of the electric DGF: 

/ii(n X Geiri, r')) = ^^ifi x Ge{r2, r')) (9a) 

n X Vi X Geiri, r') = n x V2 x Ge{r2, r') (9b) 

where ri and r2 are position vectors on either side of the boundary, /ii and ^2 are magnetic 
permeabilities on either side of the surface. Using these results, the electric field at r' can 
be expressed as a surface integral over the surface around V^. We have: 



lim 

Vs^O 



j E{r) . (V X V X Ge{r, r')) - (V x V x E{r)) . Ge(r, r') 



dr 



V-Vs 



(10a) 



jiin -^[(n5 X E{r)) . (V x Ge(r,r')) + (n^ x V x ^;(r)) . Ge(r, Ol^^^^ 



lim / [iu^^o^^{r)J^{r)-VxJ'^{r)].Ge{r,r')dr 

Vs^O J 
V-Vs 



lim 



E{r) .(hsxVx Ge(r,r') + (V x E{r)) . x Ge(r,r' 



(10b) 



dS 



Ss 



In Eq. (10) the DGF on the surface 5*5 can be split into two parts - one that is singular and 



is the DGF for infinite medium with e{r) = e{r') and ^{r) = fi{r') and the other that is 
related to the the presence of boundaries. 



Ge(r,0 = Gl°\r,0 + Gr\r,0 

=(o) =(sc) 

The equations satisfied by (r, r') and G^ (r, r') are: 

V X V X ^["\r, r') - k^V^"\r, r') = l5{r - r') 



(11) 



(12a) 



— (sc) — \ sc) 

VxVxGl {r,r')-k^Gl {r,r') 



={o) 

Gg (r, r') is given by the equation: 



where 



^?(°)(r,r') 



exp{ik\r — t'|) 
Airlr — r'\ 



(12b) 



(13) 



(14) 



At this point, the normal approach is to evaluate the surface integral in Eq. (10) by ap 



= =(o) 

proximating Ge{r,r') with G^ {r,r'). This is valid when the current densities are smooth 
functions of position[l7]. As we will see soon, the fluctuation-dissipation theory for ma- 
terials with spatially non-dispersive materials results in the current cross-spectral density 
becoming a Dirac-delta function. In such a situation, it becomes necessary to retain the 



contribution from the scattered DGF. Using Eq. (10), Eq. (12), Eq. (A5), and Eq. (A6), 
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the electric field at r' is given by: 



Eir') 



lim 



tcufxofx (r) J^^ (r) . Ge(r, r') - J"^ (r) . Gs(r, r') 



V-Vs 



lim 



1 



zuje„e ir' 



{ SC) ( SC) 



-L.J^ ir' 



dr+ (15) 



where GEir,r') = V x Ge{r, r'), and L is given by: 



L = lim (h fisV g^°\r , r') 

S6 



(16) 



The L dyadic in Eq. (15) has been discussed in great detail by Yaghjian [H]. It depends 
on the shape, and not size, of the principal volume Vs. The most important properties of 
the L dyadic that we shall make use of, without explicitly mentioning them, are that it is 
real and symmetric. The expression for the electric field is different from the expression 
found in other sources by the presence of the integral over the infinitesimal volume Vs. 
The symmetry between the equations for the electric and magnetic field suggest that the 



magnetic field must also be expressed in a same manner as Eq. (15). Indeed, that is the 



case. Just as we derived the appropriate boundary conditions for the electric DGF in the 



process of inverting Eq. (7a), so too can be done for the magnetic DGF which obeys the 



following boundary conditions: 

£i(n X G„(ri, r')) = e2{n x G^{r2, r')) (17a) 

n X Vi X G„(ri, r') = n x Va x G„(r2, r') (17b) 

Of course, the magnetic DGF, Gm(ri,r'), obeys Eq. (8). It can be seen from Eq. (9) and 



Eq. (17) that the electric and magnetic DGF are electromagnetic duals of each other |32] 



The magnetic field can then be written as: 



H(r') 



lim / Uu;eoe{r)J"'{r).GUr,r') + J''{r).GM{r,r') 

Vs^O J L 



V-Vs 



lim 

Vs^O 



1 



=(sc) (.ir) 

tcosos (r) J- (r) . G„ (r, r') + (r) . G,, (r, r') 



dr+ 



dr- 



:i8) 



iu/jioH (r') 



-L.J-(r') 



where GM{f,f') = V x Gm{r,r'). But for the L dyadic, the response functions Ge{r,r'), 
GE{r,r'), GMir.r'), and Gm{r,r') are simihar to the response functions Xee^Xeh^Xhe^ 
and Xhh defined by Agarwal [H]. The four DGFs obey the following reciprocity relations 



(obtained by using Eq. (B4)): 



-T 



/i (ra) Gg (r2, n) = /i (ri) Ge(ri, rs) (19a) 



-T 



e (ra) G„(r2, n) = e (n) G™(ri, ra) (19b) 



-T 



GMir2,n) = GEir^,r2) (19c) 



III. LDOS FROM FIELD ENERGY DENSITY AT A POINT 

To compute the energy density and the LDOS at any location, we need to compute prod- 
ucts of the type EoEpE* and fioHpH*. Since the source of this energy density is stochastic, 
the quantities of interest are cross-spectral densities (EpE*) and (HpH*), where the brack- 
ets indicate an ensemble average over all possibe configurations of the field quantities. The 
cross-spectral densities of the field quantities depend on cross-spectral densities of the elec- 
tric and magnetic current densities, i.e. terms of the form ( Jp(r'i) J|*(r'2)), {J^{fi)J^*{f2)), 
and {Jp{ri)Jjp*{r2)). The relation between the cross-spectral densities of the fluctuating 
current density and temperature can be obtained from the fluctuation-dissipation theorem, 
which states that: 

{j;{ri)r,*ir2)) = '^u6o6"Qiu,T)6{r^ - r2)(5p, (20a) 

(J™(ri)J7*(r2)) = -W'0(^,T)5(ri - r2)Spg (20b) 
{rJri)Jrir2)) = (20c) 



Using the expression for electric field in Eq. 15, an expression for eo{EpEg) in terms of 



the DGF, L, and the current densities as (terms containing ( Jp(r'i) J!"*(r2)) are neglected 



because of Eq. 20c ) 



10 



coy 
cJ 



+ Eo lim 

Vs^O 



V-Vs V-Vs 

{jr{r,)Jr{r2))GEsp{rur')G*Et,{r2,r')dr,dr2 



V-Vs V-Vs 



Vs 
Vs 



(21) 



(jj'^eols ir') r 



Substituting the results of Eq. (20a) and Eq. (20b) in Eq. (21), eo{Ep{r',uj)E*{r',uj)) can 



be written as: 



2 uj 



6o{Ep{r',uj)E;{r',u)) = --Q{uj,T)x 



Um j (5''(r)|/i(r)r5^asp(r,OGL,(r,0 +/i''(r)G^,,(r,OGL,(r',r')) dr 

V-Vs 



e"(r' 



[r 



e [v 



(22) 



+ 



LspL^q 



Unless e"{r') = 0, the last term in Eq. (22) is not finite. This can be resolved only by taking 
into account the spatial dispersion of the permittivity and permeability. We shall see how we 



can overcome this problem by using Eq. ( 11 ). The second volume integral in Eq. \\22h can be 



simplified by using Eq. (B3). Put Gi{r, r') = fi{r)Ge{r, r') and G2{r, r') = V x G^{r, r') 



GE{r,r') in Eq. (B3) to get 



lim 

Vs^O 



V-Vs 



fiir)G,ir,r')) . V x V x G,(r, r') - 



Vx (fi{r)G,{r,r'))) . VxG,(r,r' 



dr 



lim 



T 



fi{r)Ge{r,r')) . n, x V x G,(r, r') 



dS 



Ss 



(23) 
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The reason that the domain of integration of the volume integral can extend over the whole 
volume V — Vs, despite boundary surfaces in that volume of integration is because of the 
continuity conditions that the DGF satisfy (Eq. ^ for electric DGF and Eq. (17) for 



magnetic DGF). Using Eq. (11) the surface integral in Eq. (23) can be split into three 
terms as (the fourth term tends to 0): 



T 



hm^ (p {fi{r)Ge{r,r')) . [fis x V x G^{r,r'))dS 



Ss 



lim 

Ss^O 



Ss 



lim 

Ss^O 



Ss 



lim 

Ss^O 



Ss 



={o) 

/i(r)Gg (r,r' 



. (^5 X V X G^{r,r')) dS+ 



^{r)G 



r, r 



hsxV X G^ (r, r') dS+ 



={sc) 

/i(r)G^ (r,r' 



T 



.(nsxVx G„(r,r')) dS 



(24) 



The first surface integral on the RHS of Eq. (24) corresponds to the cross-spectral density 
of electric fields in the absence of all scatterers or an infinite, homogeneous medium with 
permittivity and permeability given by 6{r') and ^i{r'). This term along with the last term 



of Eq. (22) gives the cross-spectral density of electric fields in homogeneous, infinite media. 
Only when the medium is transparent is this term finite. When the medium is absorptive, 
these terms cannot be evaluated unless spatial dispersion is taken into account. In any case, 
what can always be determined is the change in the cross-spectral density because of the 



presence of scatterers. The surface integrals in Eq. (24) can be simplified using Eq. (B6) 



and Eq. (B5) to yield 



hm^ <p [^i{r)Ge{r,r')] . (n^ x V x G^{r,r')]dS 



Ss 



lim 

Ss^O 



=(o) 

/i(r)Gg (r,r' 



hs X V X G„ (r, r' 



Ss 



F*(r')= =(sc)* 
f^%r')J^L.G^ (r',r')+MO 



dS- 



=(sc) 

g: \v'y\ 



(25) 
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Equation 23 is re-written as: 



lim 

Vs^O 



T 



Mr)Ge(r,r') . V x V x G,(r, r') - 



V-Vs 



dr 



lim 



=(o) 



(26) 



dS+ 



F*(r')= =(sc)* 



= =isc) =(sc) 

L.G, {r',r')-G, (r'y) 



Taking the complex conjugate of Eq. (26) and subtracting from Eq. (26), we get: 

,2 T * T * 



lim 

Vs^O 



UJ 



5"(r) |/.(r) I^G, (r, r') . G,(r, r') + /."(r)G^(r, r') . G^(r, v') 



V-Vs 



e"(r')= ={sc)* 



er') 



e"(r')= ={sc) 



(27) 



lim (J) 3" 

Ss 



=T* I =(o) 

\x{r)G^ {r,r') .{nsxVxG^ {r,r') 



dS + Q[lx{r')G[''\r',r') 



The pq component of the LHS of Eq. 27 is the term within [ and ] in Eq. 
that the \eo{Ep{r' ,uj)E*{r' ,uj)) can be written as: 



22 



-eo{E,{r',uj)E*(r',uj)) = Q{uj,T)Ce„{r' , 



It follows 



(28) 



where Cepqir' , u) is the pq element of Ce- Ce is a cross-spectral density matrix of the electric 
fields (multiplied by £0/2 to convert units to that of spectral energy density) at the same 
point in space. Cepq{r' ,uj) is infinite in an absorbing medium when spatial dispersion of the 
permittivity and permeability are not taken into account. However, can be written as a 
sum of two terms, Ce^q{uj) (an implicit function of position, through the material properties) 
and Cepq {r',uj), an explicit function of position originating from scattering. When spatial 
dispersion is not accounted for, Cepl{uj) is finite only for a transparent medium, whereas 



Cep'^{r',uj) is finite for most cases of interest. From Eq. (27) we see that: 



=(sc) 

a (r'cu) 



TXC 



UJ I (.sr) 



and, for non-absorbing media: 



=(0) 

TXC^ 



(29) 



(30) 
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Similarly, we have for the magnetic field contribution: 

1 



2 



fio{Hp{r',u)H;{r',uj)) = Q{uj,T)Cmpg{r' ,00) (31) 
^t\r',u;) = ^3 (e{r'0^\r',r')) (32) 

TTC \ / 



and, for non-absorbing media: 

^^n!{r',u;) = —^-s (e{r'02\r',r')) (33) 

TIC \ / 

The energy density due to electric field is given by f/J = ^eo{Ep{r' ,Lj)E*{r' ,00)) = 
0(c(j, T)Tr ^Ce(r', cu) j and that due to magnetic fields is given by f/™ = 

lfXo{Hp{r',iu)H;{r',uj)) = e{co,T)Tr (^C„(r', cu) j . The trace of {r',to) and (r',cu) 
are related to the electric and magnetic LDOS. Taking the trace of Eq. (29) and Eq. (32) 
we obtain: 

C^'"^ = Tr (^^r\r', u)^ = -^Jr^ (^^(r')^rV', ^')) (34a) 

= Tr (Wt^\r\ u)^ = ^Jr^ {e{r'0^\r' , r')) (34b) 
and, for non-absorbing media: 

= Tr (^rV',^)) = ^,Tr-s (/x(r')^l°V', ^')) (35a) 

Cl?} = Tr (^l°\r',^)) = ^Jr^ (.(r')^™ (r', r')) (35b) 

To obtain the electric and magnetic LDOS in media with itself, c'f'^'^ and Cm^'' have to be 
multiplied by the appropriate factor from Eq. ^ or Eq. ([s]). 

IV. DISCUSSION 

The electric and magnetic field contribution to LDOS at a point where the medium is 
non-absorbing is given by: 

p'^ = e(r')A,TrQ ( iiir')VJr', r')] (36a) 
= p{r')^Tr^ (eir')Wmir',r')) (36b) 

TTC V / 
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The expressions for the electrical and magnetic LDOS obviously coincide with the expressions 



for LDOS in vacuum. Let us now see what the expressions for LDOS in Eq. (36) result in 
for a transparent medium of permittivity e and permeability /i. We know from statistical 
mechanics that the DOS depends on the refractive index n = J~eli as 



From Eq. (36), we have 



PsM - ^2^3 y^^i 



P = P = ^P- 



P = ^ = PsM (38) 
For a material with = 1 and e = n^, we see that the magnetic LDOS is completely from 
the magnetic field energy contribution whereas the electric LDOS has contributions from 
field energy as well as polarization energy or energy associated with the medium. 



To summarize, the concept of LDOS has been extended to material medium. A 
formal proof for LDOS in isotropic and temporally dispersive (and so absorbing) material 
is provided. In such materials, energy density, and hence the LDOS, is infinite if spatial 
dispersion is not taken into account. To include spatial dispersion in every scattering 
problem involving lossy materials makes it extremely difficult to analyze. The method 
described in this paper circumvents this problem by splitting the LDOS into a term 
that takes into account spatial dispersion by solving the problem for a homogeneous, 
infinite, spatially dispersive medium and another that does not need spatial dispersion 
but accounts for scattering from inhomogenieties, thereby simplifying the problem. It has 
been shown that the electrical and magnetic DGF that contribute to the LDOS are in 
fact electromagnetic duals of each other. This could be useful in cases where LDOS inside 
material medium become important, for instance in determining the lifetime of carriers in 
semiconductors, or to determine forces between objects surrounded absorbing media, such 
as liquids. 
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APPENDIX A: VECTOR-DYADIC RELATIONS 



A . (^B X = {A X B) . G = -B . (^A X G 

G = n ( h . G] — n X [fix G 



[F{r) . (V X V X G(r,r'))-(V x F{r)) . (V x G{r,r'))]dr 



-in. {F{r) X V X G(r, r'))dS 



(Al) 
(A2) 

(A3) 



[F{r) . (V X V X G(r, r')) - (V x V x F{r)) . G(r, Olrfr 



n . [{F{r) X V X G(r, r')) + V x F{r) x G{r, r'))]dS 



- (b[{hx F{r)) . (V X G{r,r')) + (n x V x F{r)) . G{r,r')]dS 



[F{r) . (n X V X G(r,r') ) + (V x . (n x G{r,r'))]dS 



(A4) 



jiin ^ ^;(r) .[nsxVx G^°\r, r') ] dS = -E{r') + E{r') . L 

Ss 



Proof of Eq. (A5) 



(A5) 



jiiii^ J) E{r) . Uis X V X Gt\r,r')^ dS = Xim^j E{r) . [hs x Vgo{r,r') xl)dS 

Ss Ss 



lim (p E{r) . {hsVgo{r, r') - hs . V^o(r, r')) dS 

Ss^O J 



Ss 

= E{r') .1 - E{r') 

lim / {luJ^lo^^{r)H{r)) . ( hs x G^°\r, r')\ dS = - E{r') . I + 

Ss^oJ ^ ^ ' ^ y 'J iujeoe{r') 

Ss 



(A6) 
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Proof of Eq. dAel) 



jiin ^ (F(r)) .{nsx G[°\r,r') ] dS 

Ss 



Ss- 



liin y (F(r') + (r - r') . V'F(r')) ■ (ns x ^VV^ go{r,r')dS 

Ss 

liin y ((r - r') . V'F(r')) ■ (^n^ x ^ VV^ go{r, r')dS 

Ss 

- lim ((ns.V X [{r - r') . V'F(r')] A'^^olr', r'Of^'^ 

Si-^o J 

Ss 

- jim j ns . [V X F{r')] ^Vgo{r, r')dS 



Ss 

-^[V'xF(0].I 



Put F{r) = iuj^o^i{,r)H{r) to get Eq. (A6) 



(A7) 



For Eq. (A3) and Eq. (A4) to be valid, the vector and dyad should not have singularities 



within the volume V. Refer 
taken from 



for further details. Equations (Al), (A3), and (A4) are 



APPENDIX B: DYADIC-DYADIC RELATIONS 



=(o) 

{r,r') 



47rr 



2 2 

ff I —i- h 



kr k'^r'^ 



+ 



1 + ^ - 



kr k'^r'^ 



T — —T 

nx A] .B = A . {nx B 



(Bl) 
(B2) 



[G, (r, r') . (V X G2{r, r'))-{V x Gi{r, r')f . G2{r, r')]dr = 
V (B3) 

[a^{r,r').{hxW2ir,r'))]dS 



[Gi (r, r') . (V X V X G2{r, r') - (V x V x Gi{r, r')y . G2{r, r')]dr 



[(V X Gi(r, r')f . {n x G2{r, r')) - (n x Gi(r, r')Y' . (V x G2{r, r'))]dS (B4) 



'\\T 
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Uin j> fi{r)G (r, r').{hsxV x G^{r, r') ) dS = fi{r') 

Ss 



L.G{r',r') -G{r',r' 



(B5) 



Proof of Eq. (B5) 



The proof for Eq. (B5) is similar to that of Eq. (A5) 



\im^ j> n{r)G^{r,r').i^nsxV xG {r,r')^ dS = ^*{r')—^L .G (r', r') (B6) 

Ss 



Proof of Eq. (B6) 



One way to prove this equation is to expand the scattered DGF into vector spherical waves. 



Then we can use the results of Eq. (A6) to prove Eq. (B6). Equations (B2), (B3), and (B4) 
are taken from [T2] 
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Host material 

FIG. 1: Configuration of sources. All materials are thermal sources. The host material, in general, 
is a lossless dielectric material. The scatterers are all lossy dielectric materials (this term is also 
used to refer to lossy /i). Shown in the inset is the region around r' . The volume Vs is known as 
the principal volume and Ss is the surface of V5. n is the unit normal to the surface Ss pointing 
into the volume. The volume Vs (and Ss)'\s chosen such that the permittivity and permeability at 
any point on the surface are the same as that at r' . The interfaces between any two materials are 
free of surface currents. 
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